function a = k(x1, x2, theta, l)

a = theta * exp((-1/(2*l)) * ((x1 - x2)' * (x1 - x2)));

